On the universality of the 
fluctuation-dissipation ratio in 
non-equihbrium critical dynamics 

Malte HenkeP*'^ and Gunter M. Schiitz"'^* 

'' Laboratoire de Physique des Materiaux,^ Universite Henri Poincare Nancy I, 
B.P. 239, F-545O6 Vandoeuvre-les- Nancy Cedex, France 

^Centra de Ftsica Teorica e Computacional, Universidade de Lisboa, 
Av. Prof. Gama Pinto 2, P~l 649-003 Lisboa, Portugal 

Institut fiir Festkorperforschung, Forschungszentrum Jiilich, 
D-52425 Jiilich, Germany 

February 2, 2008 



Abstract 

The two-time nonequilibrium correlation and response functions 
in ID kinetic classical spin systems with non-conserved dynamics and 
quenched to their zero-temperature critical point are studied. The 
exact solution of the kinetic Ising model with Glauber dynamics for a 
wide class of initial states allows for an explicit test of the universal- 
ity of the non-equilibrium limit fluctuation-dissipation ratio X^q. It 
is shown that the value of X^o depends on whether the initial state 
has finitely many domain walls or not and thus two distinct dynamic 
universality classes can be identified in this model. Generic ID kinetic 
spin systems with non-conserved dynamics fall into the same universal- 
ity classes as the kinetic Glauber-Ising model provided the dynamics is 
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invariant under the C-symmetry of simultaneous spin and magnetic- 
field reversal. While C-symmetry is satisfied for magnetic systems, 
it need not be for lattice gases which may therefore display hitherto 
unexplored types of non-universal kinetics. 

PACS: 05.20.-y, 05.40.-a, 05.70Jk, 64.60.Ht, 75.40.Gb 
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1 Introduction 



The understanding of the long-time behaviour of strongly interacting systems 
with many degrees of freedom and evolving far from equilibrium is an active 
topic of much current interest, see [UEllSlll! for recent reviews. Besides the 
more far-reaching aspects of disordered systems undergoing glassy behaviour, 
many of the fundamental questions of non-equilibrium statistical mechanics 
can already be studied on the conceptually simpler kinetic ferromagnetic 
systems. In several instances, general ideas can be subjected to exacting 
tests because several non-trivial and exactly soluble models are available. 

In this paper, we consider the kinetics of a purely classical spin system 
with a nonconserved order parameter and an equilibrium critical tempera- 
ture Tc > and quenched to a final temperature T < from some initial 
state. Then fluctuations of the initial state will lead on the microscopic 
level to the growth of correlated domains and the slow movement of the do- 
main boundaries will drive the irreversible time-evolution of the macroscopic 
observables. It has been realized in recent years that the associated age- 
ing effects are more fully revealed through the study of two-time correlation 
functions C{t,tw) and response functions R{t,tuj) (see section 2 for the pre- 
cise definitions) and where t is referred to as observation time while tu, is 
called the waiting time. In addition, it has been understood that correlation 
and response functions must be studied together, since out of equilibrium 
the fluctuation-dissipation theorem is no longer valid. It is usual to char- 
acterize the breaking of the fluctuation-dissipation theorem in terms of the 
fluctuation- dissipation ratio (FDR) j3] 

X(t,U:=Ti?(t,tj(^^|i^) ' (1.1) 

where at equilibrium, one would recover X{t,tw) = 1 (at T = 0, some care is 
needed to absorb T into the definition of the response function). In the scaling 
limit, where simultaneously t^ and t — t^ become large, one actually finds 
that the fluctuation-dissipation ratio X(t, t^) = X{t/t^) satisfies a simple 
scaling law in terms of the scaling variable x = t/t^, see |4 for a recent 
review. A quantity of particular interest is the limit fluctuation- dissipation 
ratio Xao defined by 

= lim f lim X(t,t^)) = lim X(x) (1.2) 
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For a quench into the disordered phase, it is known that = in generaL 
But for a critical quench with T = Tc, Godreche and Luck |j21[7| have proposed 
that Xoo should be an universal quantity. 

The evidence supporting this conjecture (which in this paper we shall 
call 'universality' for short) was based on the available results coming from 
exactly solvable kinetic spin systems quenched from a fully disordered state 
and from simulations in the 2D and 3D kinetic Ising model with Glauber 
dynamics, see [2\ for a review. Further supporting evidence in favour of the 
universality conjecture comes from field-theoretic two-loop calculations of the 
0(n)-model, again starting from a fully disordered initial state [Hlini- The 
universality of X^o has also been confirmed numerically for the 2D Glauber- 
Ising and voter models |l(Jj . 

Statements about the universahty of a physical quantity are best tested by 
varying important control parameters of suitable models and then studying 
their effects. Indeed, the role of spatially long-ranged initial correlations of 
the form 



where i/ is a control parameter, was studied in the kinetic spherical model at 
T = Tc fl] . It was shown that there exists an unexpectedly rich kinetic phase 
diagram, depending on u and the space dimension d. In most of these phases, 
either X^o = or else it is independent of z/, but in one phase Xoo was shown 
to be a function of u jTTj, thus furnishing an important qualification against 
the unrestricted universality of X^o-^ At present, it is not clear yet whether 
these results might not simply reflect a peculiarity of the spherical model. 
Therefore, we shall study here the non-equilibrium critical dynamics of ID 
ferromagnetic spin systems quenched to their critical temperature T = 0. 

In order to get analytical results, we consider in section 2 the exactly 
solvable Glauber-Ising model. We shall show that for initial correlations 
of the form ()1.3|) with u > 0, we have indeed universality of the entire 
function X{x), and thus in particular oiXoo, in agreement with the Godreche- 
Luck conjecture. However, in section 3 we study even more general initial 
states which consist of large ordered domains and then show that the scaling 
function X{x) as well as X^o are different from the ones obtained for initial 
states of the form ()1.3|) . We thereby identify two dynamic universality classes. 
These results are extended in section 4 to generalized two-state spin systems 

-'^In this phase both the space dimension d as well as the effective dimension = 2 + 1^ 
of the initial correlations are below the upper critical dimension d* . 
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evolving according to a non-conserved dynamics with detailed balance. We 
show that those which satisfy a certain global symmetry (which we call C- 
symmetry) are in one of the two universality classes found for the Glauber- 
Ising model. We also comment on the fact that in lattice gases, there is no 
need to satisfy C-symmetry which may lead to further non-universalities. We 
conclude in section 5. 

2 Glauber dynamics with correlated initial 
conditions 

We shall study the Ising chain with equilibrium Hamiltonian Ti = 
—J'^^SnSn+i where s„ = ±1 denotes the Ising spins and J > is the 
interaction strength. We consider translation-invariant initial distributions 
such that the initial magnetization mo = ( s„ ) should be different from ±1, 
i.e. there is a finite density of domain walls at initial time. In the litera- 
ture usually symmetric initial distributions with mo = are considered. In 
this section, we allow for general non-symmetric initial distributions with 
— 1 < mo < 1. The case |mo| = 1 requires separate treatment and is studied 
in section 3. In the absence of a magnetic field we assume Glauber dynamics 
for the stochastic time evolution of the spins. We set the time scale for 
individual spin flips to unity. 

2.1 Two-time correlation function 

It is convenient to write the two-time correlation function 

Cn{t + T,t) = Cn{t- r) := ( + r)so(t) ) (2.1) 

in the quantum Hamiltonian formalism, see ^3 E] for recent reviews. One 
has 

C„(t; t) = {s |<e-^°Vo^e-^°*| Pq ) (2.2) 
where a^'^'^ are the Pauli matrices acting on the n^^ site of the chain and 

Ho = l 5^(1 - <) (l - + (2.3) 

n 

is the Markov generator (stochastic Hamiltonian) for Glauber dynamics ^3] . 
Furthermore, \ Pq) is the probability vector representing the initial distribu- 
tion of spins and the constant summation vector ( s | is the left steady state. 
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The generator is constructed such that it satisfies detailed balance with re- 
spect to the equilibrium distribution at temperature T of the ID zero-field 
ferromagnetic Ising model with interaction strength J. This is achieved by 
setting 7 = tanh(2J/T) (we use units such that the Boltzmann constant 
/cb = 1)- We introduce the shorthand 

Cn := C„(0, 0) ; C„(t) := C„(t; 0) = ( s„(t)so(t) ) (2.4) 

for the initial correlations and for the equal-time correlation function, respec- 
tively. Notice that 

Cn{t) = C^nit) ; Co{t) = 1 Vt > 0. (2.5) 

For Glauber dynamics the time-evolution of the spin-expectation is linear 
and at vanishing temperature T = satisfies a lattice diffusion equation for 
a ID random walk with hopping rate 1/2 [12j. The propagator 

Gniy) = e-yiniy) (2.6) 

of the lattice diffusion equation involves the modified Bessel function /„(?/) 
[TB] and describes the probability of moving a distance of n lattice units after 
time y. Hence 

oo 
m=— oo 

which immediately yields 

oo 

Cn{t;r)= e-^In-m{r)Cm{t). (2.8) 

m=— oo 

It remains to calculate the equal-time two-point correlation function 
Cm{t)- For special initial distributions this has already been done in the 
classical paper by Glauber ^21- For our more general treatment we ob- 
serve that the total correlation function may be split into an interaction part 
C™*(t) and a correlation part C™''''(t). The latter one vanishes for uncorre- 
lated (infinite-temperature) initial states. Then 

oo 

cut) = e-2*J„(2t)+2^e-2*J|^|+,(2t) 

k=l 

oo 

+ i^M-kC^t) - hm\+k{2t)] Ck (2.9) 

fc=i 
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The two-time autocorrelation function then follows from (j2.8|) by setting 
n = 0. 

The Laplace transform of the interaction part of the two-time autocor- 
relation function has already been studied in detail by Godreche and Luck 
jH] and by Lippiello and Zannetti PZ].^ Here, we prefer to work in the more 
intuitive time-domain, which allows for an easier analysis of more general 
initial conditions. From (j2.9|) . we also define C™'(t;r) and C^°"(t;r) which 
will be calculated separately and C„(t; r) = C™^(t; r) -f C^"(t; r). Using the 
completeness property ^„ Gn = 1 of the lattice propagator we split the in- 
teraction part into a contribution which is large at early times and a second 
contribution which dominates the late-time behaviour. We find 

Cfit-, r) = e-(^+2*)/o(r + 2t) + e-Vo(r) [l - e-''lo{2t)] 

oo oo 

+4 ^ 5Z e-2*/|^l+fc(2t)e-V^(r). (2.10) 

m=l k=l 

For t, r 1 only the late-time part (containing the double sum) plays a role. 
Using the asymptotic Gaussian form 

GM = -^e-^ {1 + O {y-')) (2.11) 

of the lattice propagator we can turn the sums into integrals. Setting 

a=Jl (2,12) 



we obtain 



JO JO 

2 2 1 , , 

= 1 arctan a = — arctan—. (2-13) 

TT IT a 

The correlation function depends only on the scaling variable a. 

For the later comparison with the response function, we now make a small 
change in notation, in order to be compatible with the notations usually 



^Multispin correlators and associated response functions of the ID Glauber-Ising model 
have been studied recently in |18l IT!?] . Two-time correlators of the Ising chain in a trans- 
verse field are calculated in '20'. 
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employed jHllIIllITj. The quantity denoted 'time' t so far, we shall from now 
on call waiting time s = tw In turn, we call + r observation time and 
write t := t^ + r. We repeat the correlation function in this notation, in 
agreement with ^ [T7j 



2 / 9f 

C'^\t,t^) = -arctan./-^ (2.14) 

TT \ t - t^ 

For later use, we also record the derivative of C^^ with respect to t^i- Keeping 
t fixed, we have 

9 c^nu,, 1 l + 2a^ 1 _V2 X 1 

dt^'^' ^'''"^ 27ra(l + a2) t^ vr ^^^(a; + 1) ^^''^^ 

for either the scaling variable a or x = t/tw 

We still have to analyze the correlation part of the autocorrelation func- 
tion which has not been studied previously. In order to do so, we use the 
integral representation 

Gn{y) = ^/ cos(gn)e-^'^ (2.16) 

of the lattice propagator with the dispersion relation = (1 — cosg). This 
yields the exact expression 

oo 

CQ°"(t^;T) = e '^Im{T)C^"(tw) 

m=—OD 
oo ^ 2 

= ^^^(r)^^^- / dg sin (gm) sin (gn)e~'''*". (2.17) 

m=l n=l ^ "'-TT 

For tu],T ^ 1 the sum over m may be turned into an integral and the 
asymptotic expression (j2.1ip can be used. After some algebra we find 

Cnt^-.r) = -^EC^n fdggsin (^) iFi(i; |; aV)e-^^(^+")^ 
TT'^/^C ^^-^ Jo Vvr«>/ 

(2.18) 

where i-Fi(a; b; x) is a confluent hypergeometric series ^H]- In order to analyze 
the leading behaviour for waiting times 3> 1 we have to distinguish three 
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cases, which depend on the form of the initial correlation. We shall in this 
section consider the following form 

C^rO, 0) = C„ forn oo (2.19) 

where i/ > and B are - a priori nonuniversal - control parameters. We 
also define the unnormalized first moment 

oo 

A:=^nCn (2.20) 

n=l 

of the initial correlation function. We can now identify three distinct situa- 
tions. 



Case 1: A < oo 



Consider first the situation when the series A converges to a finite value. 
Physical examples might be antiferromagnetic alternating-sign correlations 
or for rapidly decaying ferromagnetic correlations with v > 2. For large 
waiting times tyj the leading contribution to the integral in ()2.18|) comes 
from small values of the integration variable q (which describe the long wave- 
length fluctuations of the spin system). We may then expand the sine in its 
Taylor series and the integration can be performed explicitly from the se- 
ries representation of the hypergeometric function. Summing up the infinite 
series of Gaussian integrals yields the surprisingly simple exact asymptotic 
expression, to leading order in 

Crit^^r) = --^^.^. (2.21) 

Because of the extra factor 1/t^, this is asymptotically smaller than the 
interaction part (j2.13|) . Hence the nonuniversal amplitude A does not con- 
tribute to the leading late-time asymptotics of the two-time autocorrelation 
function. 

Case 2: Slowly decaying ferromagnetic correlations 

We now consider the initial correlator (j2.19p with < z/ < 2 (such that A 
does not converge). Replacing the summation over n by an integration 
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^5^C„sin(^) ^ / dyC(yvOsing2/ (2.22) 
v^w) XytwJ Jo 

= i?t^'^/^r(l - z/) cos (^) Igr-^sign(g). (2.23) 

This yields 

(2.24) 

Also in this case the correlation part of the two-time autocorrelation function 
is asymptotically small compared to the interaction part. We conclude that 
the non-universal quantities B and which characterize the initial distribu- 
tion before the quench, do not enter into the leading late-time behaviour. 

Case 3: Partial ferromagnetic long range order 

The case i/ = corresponds to an asymptotically constant spin-spin corre- 
lation function in the initial state, mimicking (partial) ferromagnetic long 
range order. Such initial states may for example be obtained by quenching 
from a uniformly magnetized initial state to zero temperature and zero field 
and had already been studied in JT] where B = was related to the initial 
magnetization. Performing the same steps as in case 2 we find 

2B 

CQ°''{t^; t) = — arctan a. (2.25) 

TT 

This is of the same order of magnitude as the interaction contribution to the 
correlation function. For the total correlation function we obtain 

2 

Co(t^;r) = 1 - (1 - 5)- arctana. (2.26) 

71 

and therefore, for t fixed we have from ()2.15p 



^-Co(t,t.) = ^— ^— — (2.27) 



dtw TCtw 1 + X \ X — 1 

in terms of the scaling variable x = t/tw This is of the same form as in the 
uncorrelated case, but with a nonuniversal amplitude 1—B determined by the 
initial long range order. The form ()2.27|) proves universality with regard to 
details of the initial distribution of the two-time correlation function, except 
for a nonuniversal amplitude. 
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2.2 Two-time response function 

Now we consider the time evolution of the local magnetization 

5n(t) = (<(t)) = (sKe-^*|Po) (2.28) 

for an initial distribution with initial magnetization mg 7^ ±1- In zero field 
Sn{t) = mo for T = Glauber dynamics. To study the hnear response of the 
system to a small localized perturbation by a magnetic field we let an external 
field act at site of the lattice. In the quantum Hamiltonian formulation 
this perturbation of the zero-field dynamics is represented by the perturbed 
Markov generator 

H = Ho + V{h). (2.29) 

The perturbation V{h) is determined by the requirement that the full gener- 
ator H satisfies detailed balance with respect to the equilibrium distribution 



exp 



(2.30) 



of the ferromagnetic Ising system with interaction strength J and local mag- 
netic field h at site 0. This requirement, on which the usual equilibrium 
fluctuation-dissipation theorem is based, does not uniquely fix V, as differ- 
ent dynamical rules may lead to the same equilibrium distribution ()2.3()|1 . In 
jH] a heat bath prescription was used to implicitly define V. Here we follow 
more closely the philosophy of Glauber |12. and define a minimally perturbed 
dynamics by 

I L 7 . ..1 ^g-(VTK_^j^ (2.31) 



At zero temperature one has 7 = 1. We shall use the dimensionless field 
strength h/T throughout this work. 

Following standard procedures we let the field act at time and calculate 
the linear response function (in units of T) 

Rn(t,tuj) = Rn{tuj',T) = Sn(t) (2.32) 

at observation time t. As before r = t — > is the time elapsed after the 
perturbation. By expanding the full time evolution operator exp {—Ht) in 
powers of h we find from (IT^ . (ITia 

Rn{tu,-.r) = -(s|<e-^^\/'e-^*-|Po). (2.33) 
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Here V is the derivative of V with respect to h/T taken aX h = 0. Using 
(|2.7|) we see after a httle algebra that the autoresponse function {n = 0) 
factorizes 

Roit^; r) = e-Vo(r) [1 - (2.34) 

into the autopropagator Go{t) and a contribution involving the two-point 
correlation function at time tw 

To calculate the interaction part of the response function we deduce by 
analogy with (j2.1(J|) 

1 - Ci°*(t) = e-^\lo{2t) + h{2t)). (2.35) 
For large times ^ 1 we then find 

K\t^; t) = = (2.36) 

With = (x — l)/2 one obtains the autoresponse function in terms of the 
scaling variable x = t/t^. Interestingly, the same asymptotic result was 
found in jSj for heat bath dynamics. 

The calculation of the correlation part of the autoresponse function pro- 
ceeds along the same lines as the calculation of the autocorrelation function. 
We find 

Case 1: A < oo 

Here 

for large t^. Comparison with ()2.35|1 shows that we have a sub leading be- 
haviour of the correlation contribution. 

Case 2: Weakly decaying ferromagnetic correlations 

Computing the correlation as above leads to 

5 ~ ^^ ,-{u+i)/2 



cnt^) = (1 - -j t;,^'^^'^^', (2.38) 

which corresponds again to subleading behaviour. Hence initial correlations 
decaying to zero do not change the asymptotic behaviour of the autoresponse 
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function. 

Case 3: Ferromagnetic long range order 



For u = (correlations decaying to a constant value B) we obtain 



B 



(2.39) 



0rt 



w 



which is of the same order as the interaction part. Therefore 




(1 - 5) _ 1 - 5 1 



(2.40) 



7Ty/2Ttw 27rt^ a 



We see from equations ()2.27|) . ()2.40|) that the same nonuniversal amplitude 
enters the two-time correlation function and the response function respec- 
tively. 

We can now state the main result of this section: in each of the cases 



The same result was obtained in [HI for different microscopic dynamics and 
uncorrelated initial states. We note in passing that one may easily check that 
also for complete initial antiferromagnetic order Cn = (—1)" the FDR has 
the same asymptotic form. In the limit x oo we find 



in full agreement with the conjecture [HIIZI that X^o is an universal constant. 

3 Low-temperature initial states 

In the previous section we assumed a translation invariant state with a finite 
density of domain walls in the Ising system at the initial time. However, at 
very low temperatures it is more relevant to study the time evolution of an 
almost ordered system with only finitely many domain walls at the initial 



1-3 the fluctuation-dissipation ratio X = R/C does not depend on the initial 
state and is given by 





1 



(2.42) 
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time. For definiteness we consider two domain walls located at sites —L and 
L respectively of the lattice. This corresponds to the initial configuration 

i'o = -..iiiTT...Tnii... (3.1) 

where the inversions of the spins occurs at the positions —L and +L, re- 
spectively. Although these initial conditions break translation invariance, we 
have chosen the coordinate system such that refiection symmetry with re- 
spect to the origin is maintained. This is not crucial, but simplifies some of 
the exact expressions to be derived below. 

In order to calculate the two-time correlation function we use the enan- 
tiodromy relation 

= BHaB~^ (3.2) 

between Glauber dynamics and diffusion-limited pair annihilation (DLPA) 
pnil3j .^ The process DLPA describes independent random walkers hopping 
with rate 1/2 to a nearest-neighbour site and annihilating instantaneously 
upon encounter. Here Hq is the transpose of the Hamiltonian for zero- 
temperature Glauber dynamics, Ha is the Markov generator for DLPA and 
B is the factorized similarity transformation B = b^^ with the local trans- 
formation matrix 

With the enantiodromy relation (j3.2p and the initial state (j3.1|) one obtains 
in the thermodynamic limit N oo 

L 

C^At) = (s|«e-^°*|Po) = (^11 ^^e-^^V,^)- (3.4) 

k=-L 

By identifying spin up with a vacancy and spin down with a particle in 
the process of diffusion-limited annihilation the vector \m,n) is the state 
with two particles located at sites m, n and empty sites everywhere else [12] . 
Hence the calculation of the two-point correlation function is reduced to a 
two-particle problem of annihilating particles. 

The two-time autocorrelation function at site n = is given by ()2.8|) in 
terms of the equal-time correlation function. By refiection symmetry one has 

•^This relation is not to be confused with the duahty relation ^2] between Glauber 
dynamics and diffusion-limited pair annihilation. 
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C_m,o(^) = Co,mit) and therefore we may write 

oo 

Co(t^; r) = 1 - 2 ^ e-V„(r)[l - Co,^(t^)]. (3.5) 

771=1 

With (D one has 

oo oo / L \ 

i-CmW= E 5^ (s| 1- n l^'2/)^(^'i/;'^l^'^;o) (3-6) 

x=—ooy=x+l \ k=~L / 

where the two-particle propagator 

Fix, y- t\k, I- 0) = e-2* [4_,,(t)/i_,(t) - (3.7) 

for DLPA is the probabihty that two independent random walkers who 
started at fc, I have reached sites x, y after time t without having met on 
the same site. This yields 

L L+m 
y=—L x=L+l—m 

We are interested in the derivative dC{t,tw) / dtyj for large t. For an 
initial size 2L + 1 of the domain of flipped spins we consider the regime 
where ^ since at earlier times fluctuations in the initial positions of 
the domain walls have not reached the origin and hence would lead to trivial 
behaviour. To calculate the derivative of the correlator we express the exact 
expression ()3.5p in terms of the variables as in ()2.14|) . Some algebra 
similar to the previous sections yields for fixed t 

With ()2.31|) . ()2.33p the autoresponse function is given by 
R{t^-T) = Go(r)^(2-Co,i-C^-i,o(U) 

L 

= Go(r) [GL{t^) + GL+i{t^)] Gy{K) (3.10) 

y=-L 
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where we used ()3.6p . In the regime ^ this reduces to 

i?(t.;r) = ^^-i--t-3/2. (3.11) 

As expected both the autocorrelation function and the autoresponse func- 
tion contain the nonuniversal amphtude 2L + 1 which is the initial size of 
the flipped domain. However, the FDR is a universal function of the scaling 
variable x = 1 + 2a^. We find 

X(t,U =X(x) = \ (3.12) 

1 + ^/x — 1 arctan (v^; — 1) 

This scaling function is different from ()2.41|) . in particular 

X^ = (3.13) 

in contradiction to the universality hypothesis for critical dynamics as for- 
mulated in ini- 

It might be helpful to restate our results in terms of scaling functions. In 
the ageing regime, where both t^, and t — become large, one expects at 
criticality, see e.g. |2] 

Cit, s) = s-^fcit/s) , Rit, s) = s-'-'^fnit/s) (3.14) 

where a is a non-equilibrium exponent and such that for large arguments x — *■ 
OO) fc,R{x) ~ x"'^'^'^/^ which defines the autocorrelation and autoresponse 
exponents Ac and Xr, respectively. This implies that X{x) ~ x^^'^"^^-*/^ for 
X — > cxD. The dynamical exponent z = 2 throughout in the model at hand 
but the other exponents depend on the initial conditions as follows. If we 
take an initial state with decaying correlation of the power-law form (j2.22p . 
we read off from the results of section 2 

a = , \c = 1 = \r (3.15) 

However, for an initial state of the form ()3.1|) . we find 

a = ^ , Xc = , Xn = l (3.16) 

We therefore see explicitly that the different forms of the scaling function 
X{x) signal two distinct dynamical universality classes. 
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4 C-violation and nonuniversality 



In the previous sections we considered special spin flip dynamics which in 
zero field reduce to Glauber dynamics. The most general flip rates for a local 
magnetic fleld which (i) satisfy detailed balance with respect to the equilib- 
rium distribution ()2.H()j) . (ii) correspond to reflection-symmetric flnite-range 
interactions, and (iii) are nonconserved and in particular lead to Glauber 
rates for h = are as follows: 



TTT - 


-> III with rate 


(l-7)/ie-'^/^ 


Tit - 


^ III 




rn - 


- Tii 




Tii - 


- TTi 




itr - 


- iit 




lit - 


- itr 




in - 




(l + 7)/3e-^/^ 




- iti 


(l-7)/3e^/^ 



Here fi = fi{h/T) which may also depend on the neighbouring spin variables 
and must be such that /i(0) = 1. With these functions the flip rates at site 
may be written 



w{h) 



7 

[/l+2/2 + /3 + (/l 

7 



-So [S-l 



h) (s-i + ^^i) + (/i 



- 2/2 + h)s^isi + srt] 
Si) + g2S-iSi + srt] . 

(4.1) 



The short range terms SRT (which vanish for pure nearest-neighbour inter- 
actions) involve lattice sites at distances > 1 from the origin. Glauber 
dynamics [12] corresponds to Qq = 1, gi = g2 = With these rates the 
Markov generator H takes the form H = Y2n with the local spin flip 
generators hn deflned by ()2.3p for n 7^ and 



"0 = 



CTn 



)wih). 



(4.2) 



Here w{h) is the diagonal matrix obtained from 1)4.111 by replacing the clas- 
sical spin variables Si by the Pauli matrices af jTl?j . 
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To first order in h the perturbation V of tlie stocliastic time evolution 
tliat corresponds to this general choice of rates is obtained from the expansion 
of w{h) to first order in h. This leads to the following general form of the 
response function 



R{t^;r) = GoiT) 



so - 2 i^- 



+ si) [so-g'o 



9i 



-1 + si) - fi-aS-iSi - . . . 

(4.3) 

where the correlation function is evaluated at time t^. The dots denote 
contributions from the SRT next-nearest neighbour interactions and the con- 
stants g'j^ are defined by g'^ = Tdgi/dh\h=o- Glauber dynamics as used in 
the previous sections corresponds to the choice g'i = [121. The heat bath 
dynamics of corresponds to the choice g'^ = g'2 = 0, g[ = —7/2. 



4.1 C- invar iance 

In the context of spin systems the stochastic dynamics defined by the rate 
functions fi must remain invariant under simultaneous reversal of all spins 
Sk —Sk and change of sign in the external magnetic field. This global 
symmetry is an automorphism on the state space analogous to charge conju- 
gation in the quantum field theory of elementary particles. Because of this 
analogy, we shall refer to it as C-symmetry. C-symmetry of the rates requires 
the following properties 

Mh) = U-h) (4.4) 
f2{h) = U-h) (4.5) 

of the rate functions. Therefore g'oiO) = (72(0) — C-invariant systems. 

For translation-invariant initial distributions the response function (j4.3|) 
at T = (corresponding to 7 = 1) then takes the form 



1 - C,it^) + g[-^C,it^) - ( (so - I {s-i + si))(. . .) ) 

(4.6) 

where eqs. fl2.2f2.3|l were used. The dots denote further SRT terms coming 
from unspecified non-nearest-neighbour short range interactions. For long 
times the first contribution (proportional to 1 — Ci(t)) contains no nonuni- 
versal parameter of the dynamics. This is indeed the leading contribution 
since the second contribution, being a time derivative with respect to t^, is 
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subleading. Finally, the third term is a second-order (lattice) partial space 
derivative which by dynamical scaling is of the same subleading order in time 
t ime-der ivat ive . 

4.2 Consequences for non-universality 

The C-symmetry is a physical requirement that any spin dynamics must sat- 
isfy. However, the Ising model may also be regarded as a classical lattice 
gas model, see e.g. [23 1211 • ^^i^ interpretation the equilibrium distribu- 
tion ()2.3U|) describes a system of hard-core particles with attractive nearest- 
neighbour interactions which may occupy the sites of a lattice, as indicated 
by the occupation numbers nj = (1 — Sj)/2 G {0, 1}. The external magnetic 
field then corresponds to a chemical potential and a spin-flip corresponds to 
particle/vacancy interchange, i.e., an exchange of particles with an external 
reservoir. Clearly, in this interpretation of the same model there is no need 
to enforce the constraints (j4.4|) . (j4.5|) . Hence it is interesting to investigate 
universality in the absence of this symmetry. 

The response function contains the same universal part R^^^ from eq. (j4.6|) 
as in the symmetric case plus two further terms 

-Go(r) [(1 - 7)mofi'o - (7"^o - ( Si{t^)so{t^)si{t^) )) g'2] ■ (4.7) 

The first term in the second line vanishes at T = since then 7 = 1. The 
second term is non-zero only for symmetric initial states with mo = 0. As a 
function of time it is of the same order as R^'-^'^ituj; r). 

5 Conclusions 

Our study of universality in the critical, nonconserved, dynamics of purely 
classical quenched lattice models has uncovered some of the basic mecha- 
nisms behind the universality of the scaling of correlation and response func- 
tions. Explicit calculations for ID kinetic spin systems with non-conserved 
dynamics and which generalize the Glauber-Ising model have shown that the 
asymptotic form of the two-time response function does not depend on the 
microscopic form of the stochastic dynamics, provided only that C-symmetry 
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holds. ^ Furthermore, in the scahng hmit ^ 1 and t — ty^ ^ 1 both the re- 
sponse function and the time-derivative of the associated two-time correlation 
are largely independent of the form of of the correlations in the initial state, 
in agreement with the expected universality of the fluctuation-dissipation ra- 
tio X = X{x) and its limit values X^o = lim^^^oo -^(a^)- This universality 
holds true with respect to the long-range decay of the initial correlations (see 
eq. (|1.3p with u > 0) and also with respect to the 'details' of the dynam- 
ics. On the other hand, we have also shown the existence of two distinct 
kinetic universality classes which arise for initial magnetization \mo\ = 1 and 
I mo I < 1, respectively. 

In the absence of C-invariance these results only hold true for symmetric 
initial states with tuq = 0. For tuq 7^ 0, a non-universal part will contribute 
to the long-time behaviour of two-time correlation and response functions. 
Since lattice gases are in general not C-invariant, a study of these system 
will make such terms apparent. In this context, recent attempts ^31 1^ to 
construct algorithms which allow to measure the two-time response directly, 
may be of value. 

Finally, the available exactly solvable models (the ID Glauber-Ising 
model and the kinetic mean spherical model) have revealed two distinct 
possible routes towards modifications of critical dynamics beyond a fully 
disordered initial state: (1) the presence of large ordered domains and (2) 
the interplay of strong initial fluctuations with the thermal fluctuations of 
the bulk. However, in these models, the two-points functions decouple and 
can be studied independently of any longer-ranged correlations. A test of 
our scenario of possible non-universalities in more general systems, such as 
the 2D Glauber-Ising model, is called for. 
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^For the ID Glauber-Ising model, we have explicitly shown that our realization of 
Glauber dynamics reproduces the same results as the heat-bath dynamics studied in plll7|. 



20 



References 

[1] M.E. Gates and M.R. Evans (eds), Soft and Fragile Matter, Proc. 53'''^ Scot- 
tish University Summer Schools in Physics, (lOP, Bristol 2000). 

[2] C. Godreche and J.M. Luck, J. Phys. Cond. Matt. 14, 1589 (2002). 

[3] L.F. Cughandolo, |cond-mat/02103T2| 

[4] A. Crisanti and F. Ritort, J. Phys. A36, R181 (2003). 

[5] L.F. Cugliandolo, J. Kurchan, and G. Parisi, J. Physique 14, 1641 (1994). 

[6] C. Godreche and J.M. Luck, J. Phys. A33, 1151 (2000). 

[7] G. Godreche and J.M. Luck, J. Phys. A33, 9141 (2000). 

[8] P. Galabrese and A. Gambassi, Phys. Rev. E65, 066120 (2002). 

[9] P. Galabrese and A. Gambassi, Phys. Rev. B66, 212407 (2002). 
[10] F. Sastre, I. Dornic and H. Ghate, cond-mat /0308 1 78| 
[11] A. Picone and M. Henkel, J. Phys. A35, 5575 (2002). 
[12] R. Glauber, J. Math. Phys. 4, 294 (1963). 

[13] G.M. Schiitz, Exactly solvable models for many-body systems far from equilib- 
rium, in: Phase Transitions and Critical Phenomena, Vol. 19 Eds. C. Domb 
and J. Lebowitz (Academic Press, London 2000). 

[14] M. Henkel, in A. Kundu (ed), Classical and Quantum Nonlinear In- 
tegrable Systems: Theory and Application, (lOP, Bristol 2003), ch. 10 
(co nd-mat/030351 2). 

[15] B.U. Felderhof and M. Suzuki, Physica 56, 43 (1971). 

[16] I S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Series and Products 
(Academic Press, Orlando 1981). 

[17] E. Lippiello and M. Zannetti, Phys. Rev. E61, 3369 (2000). 

[18] P. Mayer, L. Berthier, J.P. Garrahan and P. Sohich, Phys. Rev. E68, 016116 
(2003). 

[19] P. Mayer and P. Solhch, |cond-niat/03d72T4i 



21 



[20] V. Subrahmanyam, |cond-n iat/0308393l 

[21] M. Henkel, E. Orlandini and G. M. Schiitz, J. Phys. A28, 6335 (1995). 
[22] J.E. Santos, J. Phys. A30, 3249 (1997). 

[23] R. Kubo, Statistical Mechanics (North Hohand, Amsterdam 1965). 

[24] B. Chopard and M. Droz, Cellular automata modelling of physical systems 
(Cambridge University Press, Cambridge 1998). 

[25] C. Chatelain, [conci-mat/0303545[ 

[26] F. Ricci-Tersenghi, icond-mat/03075651 



22 



